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Abstract.  Two  “non-rotating  origins”  were  defined  by  the  lAU  in  2000 
for  the  measurement  of  Earth  rotation:  the  Celestial  Ephemeris  Origin 
(CEO)  in  the  ICRS  and  the  Terrestrial  Ephemeris  Origin  (TEO)  in  the 
ITRS.  Universal  Time  (UTl)  is  now  defined  by  an  expression  based  on 
the  angle  9  between  the  CEO  and  TEO.  Many  previous  papers,  e.g., 
Capitaine,  Guinot,  &  McCarthy  (2000),  developed  the  position  of  the 
CEO  in  terms  of  a  quantity  s,  the  difference  between  two  arcs  on  the 
celestial  sphere.  A  similar  quantity  s'  was  defined  for  the  TEO. 

As  an  alternative,  a  simple  vector  differential  equation  for  the  po¬ 
sition  of  a  non-rotating  origin  on  its  reference  sphere  is  presented  here. 
The  equation  can  be  easily  numerically  integrated  to  high  precision.  This 
scheme  directly  yields  the  unit  vector  of  the  CEO  in  the  ICRS,  or  that  of 
the  TEO  in  the  ITRS,  as  a  function  of  time.  This  simplifies  the  deriva¬ 
tion  of  the  main  transformation  matrix  between  the  ITRE  and  the  ICRS. 
The  directness  of  the  development  may  have  pedagogical  and  practical 
advantages  for  the  vast  majority  of  astronomers  who  are  unfamiliar  with 
the  history  of  this  topic. 


lAU  resolution  B1.8  from  the  2000  General  Assembly,  which  establishes  the 
Celestial  Ephemeris  Origin  (CEO)  and  the  Terrestrial  Ephemeris  Origin  (TEO), 
and  redefines  UTl  in  terms  of  the  angle  9  between  them,  has  not  received  wide 
publicity  within  the  general  astronomical  community.  Among  those  who  are 
aware  of  it,  probably  few  appreciate  its  basis  or  implications.  The  CEO  and 
TEO  are  specific  examples  of  the  concept  of  a  “non-rotating  origin”  (NRO), 
first  described  by  Guinot  (1979).  The  concept  itself  is  quite  simple,  but  the 
mathematical  details  of  its  practical  implementation,  as  usually  presented  (e.g., 
Capitaine  et  al.  2000),  may  discourage  non-specialists.  This  paper  describes 
an  alternative  mathematical  development,  which  is  precise  enough  for  practical 
computations,  but  avoids  some  of  the  untidy  aspects  of  an  analytical  approach 
based  on  spherical  trigonometry.  In  the  space  allotted  here,  only  an  outline  of 
the  scheme  can  be  presented. 

A  non-rotating  origin  can  be  most  simply  described  as  a  point  on  the  moving 
equator  whose  instantaneous  motion  is  always  orthogonal  to  the  equator.  If  n(t) 
is  a  unit  vector  toward  the  instantaneous  pole  and  x(t)  is  a  unit  vector  toward 
an  instantaneous  non-rotating  origin,  a  simple  geometric  construction  based  on 
this  definition  yields  the  following  differential  equation: 


±{t)  =  -(x(t)  •  n(t))  n(t) 
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That  is,  if  we  have  a  model  for  the  motion  of  the  pole,  n(t),  the  path  of  the 
non-rotating  origin  is  described  by  x(t),  once  an  initial  point  on  the  equator, 
x(to),  is  chosen.  Conceptually  and  practically,  it  is  simple  to  integrate  this 
equation,  using,  for  example,  a  standard  4th-order  Runge-Kutta  integrator.  For 
the  motions  of  the  real  Earth,  fixed  step  sizes  of  order  0.5  day  work  quite  well, 
and  the  integration  is  quite  robust.  Since  this  is  really  a  one-dimensional  problem 
carried  out  in  three  dimensions,  two  constraints  can  be  applied  at  each  step: 
|x|  =  1  and  X  •  n  =  0. 

By  numerically  integrating  this  equation,  we  obtain  a  time  series  of  unit 
vectors,  x(t).  Depending  on  our  model  of  the  motion  of  the  pole,  n(t),  these 
vectors  point  toward  either  the  celestial  or  terrestrial  non-rotating  origin.  That 
is,  the  vectors  define  the  directions  of  either  the  CEO  in  the  ICRS  or  the  TEO 
in  the  ITRS.  If  we  are  working  in  the  celestial  frame  and  the  CEO  is  the  non¬ 
rotating  origin  of  interest,  the  model  of  the  pole’s  motion  is: 

n{t)  =  F  P"(t)  N"(t)  0  j  (2) 

where  F,  P,  and  N  are  the  standard  matrices  for  frame  bias,  precession,  and 
nutation,  respectively.  The  transpose  symbols  on  the  precession  and  nutation 
matrices  indicate  that  they  transform  vectors  from  time  t  to  J2000.0.  The  frame 
bias  matrix  transforms  vectors  from  the  dynamical  system  to  the  ICRS. 

Once  the  integration  is  completed,  the  main  part  of  the  ITRS-to-ICRS 
transformation  (which  converts  a  terrestrial  vector  to  the  equivalent  celestial 
vector)  can  be  simply  expressed,  for  any  epoch  ti,  in  terms  of  the  components 
of  three  unit  vectors:  x(ti),  n(ti),  and  y{ti)  =  n(tj)  x  x(tj).  These  three  vec¬ 
tors  have  components  expressed  with  respect  to  the  ICRS  axes  and  define  the 
orthonormal  basis  for  what  has  been  called  the  “intermediate”  coordinate  sys¬ 
tem.  This  system  has  the  celestial  pole  in  the  z-direction  and  the  CEO  in  the 
x-direction,  with  the  instantaneous  equator  as  the  xy-plane.  The  transformation 
matrix  that  takes  a  vector  from  the  intermediate  system  to  the  ICRS  is  simply: 


But  to  apply  this  matrix,  we  must  first  get  the  terrestrial  vector  into  the  in¬ 
termediate  system.  A  vector  in  the  ITRS,  after  correction  for  polar  motion,  is 
transformed  to  the  intermediate  system  by  a  simple  rotation  through  the  angle 
9,  which  is  a  linear  function  of  UTl  (the  expression  is  specified  in  lAU  resolution 
B1.8  of  2000).  So  the  complete  transformation  from  the  ITRS  to  the  ICRS  is 

r,  =  C  Ksi-e)  W'  rt  (4) 


where  is  a  vector  in  the  ITRS,  Vc  is  the  corresponding  vector  in  the  ICRS, 
W'  is  the  polar  motion  matrix,  and  R3  represents  a  simple  rotation  about  the 
z-axis.  (See  Chapter  5  of  the  lERS  Conventions  (2003)  for  the  expression  for  0 
and  the  form  of  W',  there  referred  to  as  W(t).) 

This  approach  also  provides  a  very  simple  equation  for  apparent  sidereal 
time.  Greenwich  apparent  sidereal  time  (CAST)  is  simply  the  Greenwich  hour 
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Figure  1.  Differences  in  results  of  ITRS-to-ICRS  transformation  methods: 
CEO-based  minus  equinox-based.  A  numerical  integration  of  the  CEO  posi¬ 
tion  in  the  ICRS  was  used  for  the  CEO-based  method. 


angle  (GHA)  of  the  true  equinox  of  date.  Using  the  form  of  eq.  (2),  the  direction 
of  the  true  equinox  in  the  ICRS  is  given  by 


r(t)  =  F  p^(t)  N^(t) 


(5) 


and  its  direction  in  the  intermediate  system  must  be  ( T  •  x,  T  •  y,  0 ).  If  we 
define  the  Greenwich  meridian  to  be  the  plane  that  passes  through  the  TEO  and 
the  poles,  then  Greenwich  apparent  sidereal  time  is  just 


CAST  =  GHA  T 


9  —  arctan 


X 


(6) 


since  6  is  the  angle  between  the  TEO  and  the  CEO,  and  the  latter  defines  the 
direction  x. 

Does  this  whole  approach  work  in  practice?  If  so,  its  results  must  match 
those  given  by  a  conventional  equinox-based  development.  Eigure  1  shows  the 
differences  in  the  results  of  CEO-based  and  equinox-based  methods  of  perform¬ 
ing  the  ITRS-to-ICRS  transformation  on  a  vector  in  the  equatorial  plane.  Eor 
these  computations,  the  numerical  integration  of  the  CEO  position  was  based 
on  the  motion  of  the  pole  defined  by  the  lAU  2000A  precession-nutation  model, 
as  given  in  the  lERS  Conventions  (2003).  The  new  lERS  formula  for  sidereal 
time  (including  all  the  “complementary  terms”)  was  used  in  the  equinox-based 
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Figure  2.  Same  as  Fig.  1,  except  that  lERS  algorithms  based  on  analytic 
series  were  used  for  the  CEO-based  method. 


(conventional)  transformation.  The  initial  CEO  position  used  in  the  integration 
was  adjusted  so  that  the  difference  between  the  two  methods  was  zero  at  2003.0. 
(The  ICRS  right  ascension  of  the  CEO  at  J2000.0  turned  out  to  be  0^002012.) 
Polar  motion  and  AT  were  assumed  zero  throughout.  The  figure  shows  that 
the  two  transformation  schemes  match  to  within  a  few  microarcseconds  from 
1800-2200,  with  some  divergence  outside  of  that  time  frame.  A  similar  com¬ 
parison  of  the  two  kinds  of  sidereal  time  formulas  yields  differences  that  are 
identical  by  eye  to  the  top  plot  of  Eig.  1.  The  same  systematic  pattern  (with 
some  additional  numerical  noise)  is  seen  in  Eig.  2,  which  shows  the  same  kind 
of  ITRS-to-ICRS  transformation  comparison,  but  using  only  lERS  algorithms 
throughout.  Therefore,  the  CEO  integration  approach  described  here  seems  to 
be  at  least  as  precise  as  that  provided  by  the  lERS. 

The  CEO  integration  development  and  results  will  be  described  in  more 
detail  in  Kaplan  (2004). 


References 

Capitaine,  N.,  Guinot,  B.,  &:  McCarthy,  D.  D.  2000,  A&A,  355,  398 
Guinot,  B.  1979,  in  lAU  Symp.  82,  Time  and  the  Earth’s  Rotation,  ed.  D.  D. 

McCarthy  &  J.  D.  H.  Pilkington  (Dordrecht:  Reidel),  7 
lERS  Conventions  (2003),  http://maia.usno.navy.mil/conv2000.html 
Kaplan,  G.  2004,  in  preparation 


